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ABSTRACT: We make an complementary investigation of the primordial trispectrum from exchanging intermediate scalar 
modes in multi-field inflationary models with generalized kinetic terms. Together with the calculation of irreducible con- 



tributions to the primordial trispectrum in Ref.[104|, we give the full leading-order primordial trispectrum in generalized 
multi-field models. 
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1. Introduction 



One of the most exciting ideas of modern cosmology is inflation [|J, which can solve the flatness, the horizon, and the 
monopole problem of the standard big bang cosmology. Such a period of cosmological inflation can be attained if the 
energy density of the universe is dominated by the vacuum energy density associated with the potential of some scalar 
field(s). Over the years, inflation has become so popular because of its prediction of nearly scale-invariant primordial 
density perturbation. In the inflationary scenario, the primordial fluctuations of quantum origin were generated and frozen 
to seed wrinkles in the Cosmic Microwave Background(CMB) [§] ^ [£| [|| [|] [0] [[| [§] and today's Large-scale Structure 
(LSS) (15 {J |@ [01 0. 

Inflation is mostly a framework of theories rather than a single model or theory. From the observational point of 
view, many inflationary models are "degenerate". Measuring tensor modes in the CMB anisotropy and the spectral index 
of the power spectrum of adiabatic perturbation are not adequate to efficiently discriminate among different inflationary 
scenarios. Fortunately, we have another observable available, which proves to be valuable in providing us with additional 
information beyond the power spectrum to discriminate models. It is the deviation from a purely Gaussian statistics among 
CMB anisotropies [15| Jl6[|, which arises from interaction(s) among perturbations, leading to non-vanishing higher-order 
correlated functions. Due to its importance, constraining and predicting primordial non-Gaussianity has become one of 
the major efforts in modern cosmological community. 

The simplest single-field slow-roll inflation models, within the context of Einstein gravity and the standard initial 
adiabatic vacuum, is only able to generate negligible amount of non-Gaussianity Jutl, which is undetectable by current 
observations of the CMB or even LSS. In the theoretical aspect, there are several ways to approach large non-Gaussianity. 
A short list of these models and mechanisms includes fc-inflation or models with general non-canonical kinetic terms 
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Since much more observational data will be available in the near future from WMAP/PLANCK and LSS experiments, 
it is very necessary to study the four and higher-point correlation functions. In this paper, we make a complement to the 
calculation of Ref. [104|, in which we calculated the contributions to the primordial trispectrum in general multi-field 
inflation from the irreducible or so-called "contact" diagrams. A complete calculation of the trispectrum should also 
include the contributions from reducible or so-called "exchanging intermediate scalar modes" diagrams, as performed in 
[ p2| , 105, [30| in the investigation of the trispectrum in single-field and multi-field inflationary models, and in [ |28| | 
where exchanging gravitons was considered. In this paper we show that, the contributions to the final trispectrum arising 
from exchanging scalar modes has the same magnitude as those from the contact contributions, and thus is also very 
important. 

The remainder of this paper is organized as follows. In Sec. 2, we briefly review the background evolution and linear 
perturbations for our model. Readers who are interested in the details are encouraged to refer to [104]. In Sec. 3, we 
calculate the tri-spectrum which originating from correlating (or exchanging) scalar modes. The full trispectrum, which 
includes both contacting and correlating scalar contributions, is also discussed. 



2. Basic Setup 



2.1 Model and Background 

In this work we consider a general class of multi-field models containing J\f scalar fields coupled to Einstein gravity. The 
action takes the form 



S = d 4 x 



-9 



where (fi 1 (I = 1, 2, • • • , AO are scalar fields acting as inflaton fields, and 



X 



i j 



(2.1) 



(2.2) 



is the kinetic term (matrix), is the spacetime metric tensor with signature (— , +. +, +). J"-indices are raised, 
lowered and contracted by the A^-dimensional field-space metric Gu = Gu{4> 1 ). This form of the Lagrangian in- 
cludes multi-field k-inflation and multi-DBI models as special cases. For example, multi-field k-inflation has the scalar- 
field Lagrangian as P(A,0 / ), where X = tvX IJ = GuX IJ , while in multi-field DBI models, P(X IJ , (j) 1 ) = 
-T^ry (VV - l) - V(^) with V = 1 - 2fGuX IJ + 4fX [ /xj ] - 8f 3 X [ / XjX% ] + 16f*X [ /xjX« X?. 



We work in the ADM formalism of gravitation, in which the spacetime metric is written as 

ds 2 = -N 2 dt 2 + h ij (dx i + N l dt){dx 3 + N j dt) , 



(2.3) 



where N — N(t, x) is the lapse function, iVj = Ni(t, x) is the shift vector, and hij is the spatial metric on constant time 
hypersurfaces. The ADM formalism is convenient because the equations of motion for N and N l are exactly the energy 
and momentum constraints which are easy to solve. Under the ADM formalism, the action (2.1) can be written as (up to 
total derivative terms) 



S = J dtd 3 xVhN Q# (3) + ^ ( E ij Eii ~ E2 )^j+ J dtd 3 xVhNP, 



where h = det hij and the symmetric tensor 



Ej 



(2.4) 



(2.5) 



with Vj the spatial covariant derivative defined with the spatial metric hij and E = trEij = h^Eij. is the three- 
dimensional Ricci scalar which is computed from the spatial metric hij. In the ADM formalism, spatial indices are raised 
and lowered using hij and h % i . 

In the ADM formalism, the kinetic matrix X IJ can be written as 



X 



ij 



2N 2 



(2.6) 



where v 1 = 6 1 - N^jch 1 . 
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2.1.1 Equations of Motion 

The equations of motion for the scalar fields are 

V fl (P {IJ) d^ I )+Pj = 0, (2.7) 
where V M is the four-dimensional covariant derivative. Here and in what follows, we denote 

dP d 2 P 

as a shorthand notation. 

The equations of motion for N and Ni are the Hamiltonian and momentum constraints respectively, 



+2P- Ap^V - -L (EijEV -E 2 )=Q, 



N 2 >\* J > N 2 

P 

~N 



(2.9) 



2.1.2 Background 



In this work, we investigate scalar perturbations around a flat FRW background, the background spacetime metric takes 
the form 

ds 2 = -dt 2 + a^t^cfccW , (2.10) 
where a(t) is the so-called scale-factor. The Friedmann equation and the continuity equation are 

H 2 ^^ = -(2X IJ P (IJ) -P) , 

3 3 V ' {U) ' ' (2.11) 

p = -3H(p + P). 

In the above equations, all quantities are background values. From the above two equations we can also get another 
convenient equation 

vIJ r. 

The background equations of motion for the scalar fields are 



H = -X IJ P UJ) . (2.12) 



P,{U)¥ + (3HP {IJ) + P (/J) ) ^ - P.j = , (2.13) 

where P j denotes derivative of P with respect to (j) 1 : P j = 

In this work, we investigate cosmological perturbations during an exponential inflationary period. Thus, from (2.12) 
it is convenient to define a slow-roll parameter 

e = -4 = ^#°. (2.14) 
H 2 2H 2 V ' 

2.2 Perturbation Theory in the Spatially-flat Gauge 



The scalar metric fluctuations about our background can be written as (see [106 107] for nice reviews of the theory of 
cosmological perturbations) 

SN = a, 

5Ni = di!3, (2.15) 
Sg zj = -2a 2 (ipSij - d.djE) 

where a, (3, tp and E are functions of space and time 1 . The scalar field perturbations are denoted by 5(f) 1 = Q 1 . 

Before proceeding, we would like to analyze the (scalar) dynamical degrees of freedom in our system. In the begin- 
ning we have Af + 4 apparent scalar degrees of freedom. The diffeomorphism of Einstein gravity eliminates two of them 2 , 



'This form of ansatz corresponds to <5goo = 1 — N 2 + NiN' and <5<?oi = Ni. 
2 See rtlOg] for a detailed discussion on the gauge issue of cosmological perturbations. 
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leaving us TV + 2 scalar degrees of freedom. Furthermore, two of these TV + 2 degrees of freedom are non-dynamical. 
In the ADM formalism, these are just the fluctuations 8N = a and 5Ni = di(3. Thus, there are TV propagating degrees 
of freedom in our system. As has been addressed, the diffeomorphism invariance allows us to choose convenient gauges 
to eliminate two degrees of freedom. In single-field models, there are two convenient gauge choices: comoving gauge 
corresponding to choosing 5<fi = E = or spatially-flat gauge corresponding to %jj = E = 0. In the multi-field case, the 
comoving gauge loses its convenience since we cannot set Sp = for every field in multi-field case. Thus, in this work 
we use the spatially-flat gauge. 

In the spatially-flat gauge, propagating degrees of freedom for scalar perturbations are the inflaton field perturbations 
Q 1 (t, x), while SN and <57V^ are non-dynamical constraints. In this work, we focus on scalar perturbations. In general, it 
is well-known that in the higher-order perturbation theories, scalar/vector/tensor perturbation modes are coupled together. 
However, from the point of view of the perturbation action approach, these couplings are equivalent to exchanging var- 
ious modes. In this work, we focus on interactions of scalar modes themselves, and neglect tensor perturbations. The 
perturbations take the form 

<j> I (t,x) = I o (t) + Q I (t,x), 

N = l + ai +a 2 -\ , 

N i = d i (fi 1 +fo + '-) + 6 li +6 2 i + -~ , 

where 4>l(t) is the background value, and a n , /3 n , 9 n i are of order 0(Q n ). 

The next step is to solve the constraints a n , j3 n and 9 n i in terms of Q 1 . Fortunately, in order to expand the action to 
third-order in Q 1 , the solutions for the constraints up to the first-order are adequate. At the first-order in Q 1 , a particular 



solution for equations (2.9) is: 

1 



a 



—d~ 2 

2H 



(P {IJ) +2X KL P {IJ){KL) ) (^-P^KQt-tf&j ^ 
-3HP AIJ) ^Q J - P AU)K Q K 2X IJ + PjQ 1 ' , 



hi = . 



Here and in what follows, repeated lower indices are contracted using Sij, and d 2 = didi. d 2 is a formal notation and 
should be understood in fourier space. 

2.3 Linear Perturbations 

In multi-field model, we can decompose the perturbation into one instantaneous adiabatic sector and one instantaneous 
entropy sector. The "adiabatic direction" corresponds to the direction of the "background inflaton velocity" 

e{= , V = = g-, (2.18) 



where we define a = y P^jk)4> j 4> K ■> which is the generalization of the background inflaton velocity. Actually a is 
essentially a shorthand notation and has nothing to do with any concrete field. Note that a is related to the slow-roll 
parameter e as a 2 — 2H 2 e. 

We introduce (TV — 1) basis e^, (n = 2, ■ ■ ■ ,TV) which are orthogonal with e{ and also with each other. The 
orthogonal condition can be defined as 

P.iU^t^Smn- (2.19) 

Thus the scalar-field perturbation Q 1 can be decomposed into instantaneous adiabatic/entropy basis: 

Q 1 = e^Q m , m=l,---Af. (2.20) 
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Up to now our discussion is rather general, without further restriction on the structure of P(X IJ , 
we consider a general class of two-field models, with the following Lagrangian of the scalar fields 3 : 



P(X ZJ ,tf) = P(X,Y,< 



b 1 ). In this work, 



(2.21) 



with X = Xf — GijX ij and Y = XjXf. This form of Lagrangian not only is the most general Lagrangian for 
two-field models and thus deserves detailed investigations, but also can make our discussions on the non-Gaussianities in 
two-field models in a more general background. 

After performing the decomposition into instantaneous adiabatic/entropy modes, at the leading-order, the second- 
order action for the perturbations takes the form 4 



with 



S. 



(main) 



r 3 3 /i . . i \ 

/ dtd X d i ~^K"mnQmQn ~ ^ 2 ^ rnn ^Q m ^Q n J ' 



( 



p 



(MAT)' 



iM IN 



n I K J L 

P(IK)(JL) e l e n e l e m : 



1 



SlmSln) , 



(2.22) 



(2.23) 



where we introduce 5 



P 



x 



c 2 = 



Px + 2X (P xx 
P.x 



4XP iXY + 3P r + AX^P^yy) 



(2.24) 



P x + 2XP tY 

which are the propagation speeds of adiabatic and entropy perturbations respectively. It is useful to note that K mn is 
diagonal, K\\ = l/c 2 , K.22 = l/ c e anc ^ ^12 = ^-21 = 0, as a consequence of the adiabatic/entropy decomposition. 
c a 7^ c e is a generic feature in multi-field models; this can be seen explicitly from the definitions in (2.24), the speed of 



sound for the adiabatic mode and the entropy mode(s) have different dependence on the P-derivatives 6 . 
At this point, it is convenient to introduce two parameters: 

c _ X(P XX + 2P XY ) 



P 



x 



2XP tY 
2 



A = X A P XX + -X A P. XXX + 2 i^YPy + 6Y'P. YY + — Y Pyyy 
+ 4 (X 2 YP XXY + 2XYP XY + 2XY 2 P XYY ) , 

where all quantities are background values, and we have used Y = X 2 . As we will see later, although the ^in- 
dependences of P(X, Y, (f) 1 ) in general can be complicated, the non-linear structures of P affect the trispectra through the 
above specific combinations of derivatives of P. 

After introducing new variables whose kinetic terms are canonically normalized 



a 

— 1 

c 



Qs — Qs ) 

C e 



and changing into comoving time defined by dt = adrj, the quadratic action takes the form 



5' 2 



drjd x ■ 



1 



Q' 2 + {H 2 + W) Ql - cl(dQ a f + Q' 2 + (H 2 + H') Q 2 S - c 2 (dQ s 



(2.25) 



(2.26) 



3 This form of Lagrangian is motivated from that, for multi-field fc-inflation models ^fj], the LagTangian is simply P(X, (j) 1 ). In [ ft^l a special 
form of the Lagrangian P(Y, 4, 1 ) with Y = X + (x 2 - XuX IJ ) was chosen in the investiga tion o f bispectrua in two-field models, which 

is motiv ated by the multi-field DBI action. In this work, we use the more general form of the Lagrangian ( 2.21 ). 

4 In (2.22) we neglect the mass-square terms as MmnQ m Q n and the friction terms such as ~ QmQn- In general these terms may become 
important, especially they may cause non-vanishing cross-correlations between adiabatic mode and entropy mode around horizon-crossing. See | p3| | for 
detailed investigation of these cross-correlations for the same model in this paper, and [ p6[ |s7| for recent studies on multi-field perturbations. 

5 We use c a and c e rather than c CT and c s in order to avoid possible confusion, since in the literatures c s has special meaning, i.e. the speed of sound 
of perturbation in single-field models. 

6 This fact was first point out apparently in [ p9[ |6o| ] in the investigation of brane inflation models. See also |o"i| , |6^ . [tlj for extensive 

investigations on general multi-field models with different c a and c e . 
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The action (2.22) or (2.26) describes a free theory. Performing a canonical quantization, we write 

Q a (k,ri) = akUk(v) +aL k ul(r)) , Q s (k,r]) = a k v k (i]) + a}_ k v* k {ri) , (2.27) 

where Uk {if) and v k (rf) are the mode functions, which satisfy the corresponding classical equations of motion 

u'l + [c 2 ,k 2 - (H 2 + %')] u k = , v'l + [elk 2 - {H 2 + %')] v k = . (2.28) 

Finally, what we are interested in are the tree-level two-point functions for Q„ and Q s , defined as 



(Qa(ki,m)Qa(k 2 ,m)) = (2n) 3 6 2 (k 1 + k 2 )G kl (771,772) 
(Q s (k 1 ,m)Q s (k 2 ,r 1 2)) = (27rfS 2 (k 1 + k 2 )F kl (771,772) 



(2.29) 



with 

Gk(vi, m) = u k(m) u Uv2) , F k (T] 1 ,V2) = v k (m)vt(v2), (2.30) 

where Uk (77) and Vk (77) are the mode functions for adiabatic perturbation and entropy perturbation respectively: 



Vk ( v )^^=(l + iCe kr 1 )e-^ k \ 



(2.31) 



The so-called "power spectra" for adiabatic and entropy perturbations are defined as P a (k) = Gk(i]*,f]*) and 
P s (k) = Fk{f}*, 77*), where 77* can be chosen as the time when the modes cross the sound-horizon, i.e. at c d k = aH 
for adiabatic mode and c e k = aH for entropy mode(s) 7 . In the so-called comoving gauge, the perturbation Q a is di- 
rectly related to the three-dimensional curvature of constant time space-like slices. This gives the gauge-invariant quantity 
referred to as the "comoving curvature perturbation": 

n = ^-Qa ■ (2.32) 
a 

The entropy perturbation Q s is automatically gauge-invariant by construction. It is also convenient to introduce a renor- 
malized "isocurvature perturbation" defined as 

S=—Q s . (2.33) 
a 

In the cosmological context, it is also convenient to define the dimensionless power spectra for comoving curvature 
perturbation and isocurvature perturbation respectively: 

H 2 „ H 2 k 3 „ 1 /m 2 



a to 2ee a \M 

H 2 „ H 2 k 3 „ 1 / H^ 2 



Vs * ~ & 2 Vs * ~ a 2 27r2 Ps * (fc) 2ec e 
In the above results, all quantities are evaluated around the sound-horizon crossing 



3. Non-linear perturbations 

In this section, We calculate the tri-spectrum which comes from correlating (or exchanging) scalar modes. The full 
trispectrum which includes both contacting and correlating scalar contributions is also discussed in this section. 

7 In general multi-field models, adiabatic/entropic modess with the same comoving wavenumber k exit their sound-horizons at different time, due 
to their different speeds of sound, c a 7^ c c . This fact will bring new interesting phenomenology in multi-field models. As was shown in the 
cross-correlations between adiabatic/entropic modes would be enhanced by a small c e / c a ratio. 
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3.1 Trispectra from Correlating Scalar Mode 



The third-order action for the model (2.1) has been derived in [WM: 



s: 



(main) 



dtd X 0, [ ^ m n I Q m Q n Q I n n T mnl Qm^iQn^iQt 

2 2a z 



(3.1) 



with 



i — 'mnl 



P,{IK){JL)Z\Z m e n e l + o ( P ,{MN) 4^ 4" ) P ,{IK){JL){PQ) 



1M IN 



,p I p K p J p L p P p Q 

\e 1 e m e 1 e„e 1 e x 



(3.2) 



P(MN)<f> M <t> N P(IK)(JL) e[e*e J n ef 



In this article, we still work on the double-field model. It is a straightforward task to generalize our calculation to a more 
general multi field model. 

Direct algebra gives the cubic-order interaction Hamiltonian: 



Hi(t) = / drd 3 x 



~Z a Q'* + ^T„Q' a d i Q«d i Q 



a 



^cQ'M's + I'T.C&diQAQs + ^cQsdiQadiQ. 



(3.3) 



where the subscript "I" denotes the interactional picture and the five effective couplings S CT etc. are given in Appendix A. 



The trispectrum is the four-point correlation function of perturbations. According to the in-in formalism [ 108 1, The 
trispectrum which comes from scalar exchanging can be formulated as 



(Q 4 ) D -23ft 



drf / dn"(0 I \QjH I (r,')H I (r 1 ")\0 I ) 



drf \ dr 1 "(0 I \H I ( V / )QjH I (r 1 ")\0 I ) 



(3.4) 



The calculation is straightforward but rather tedious. Here we simply collect the final results. The leading contribution 
from exchanging an intermediate scalar mode to the purely adiabatic four-point function (Q 4 -) is given by (see Appendix 
H for details) 



(Qa (T, fcl) Qa (T, k 2 ) Qa (T, fe 3 ) Qa (t, fe^SE 

4 -9. 



= (27r) 3 5 3 k i ) { 2 S ^ ^a^a (c a kl,Cak2,Cak 3 ,Cak4,C a k 12 ) 

- 3E a T^c 8 a 1^ (c a k ll c a k 2 ,c a k 3 ,c a k 4 ,c a k 12 ) + 22 { c 2) {c a k\, c a k 2 , c a k 3 , c fc 4 , c fei 2 ) + 23perms| , 

(3.5) 

where "23 perms" denotes the other 23 permutations among four external momenta fei, • • ■ , k 4 . In (3^), the integrals I a , 



-^t^ ( c aki,c a k 2 ,c a k 3 ,c a k4,Cak 12 ) (c a k 1 ,c a k 2 ,c a k 3 ,c a k4,c a k 12 ) + 1^ (c a ki,c a k 2 ,c a h,c a ki,c a ki2) 
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etc are defined in Appendix ^. The mixed adiabatic/entropy four-point function (Q^Qs) i s gi ven by 



(Q<t (r, fei) Q CT (t, fc 2 ) Q s (r, fe 3 ) Q s (r, A;4)) SE + 5perms 



1 

:(2n) 3 S 3 (^2 fe «){ 3S o- S cCaC*I (c a fci,C Q fc 2 ,C e fc3,Cefc4,C Q fci2) + 2S^C^C^I a (c a fcl,C e fc3,C a fc 2 ,Cefe4,C e fcl3) 

z 

+ T CT T s c^Cg p b (1) (c a fci,c Q fc2,c e fc 3 ,c e fc4,c a fci2) +2I^ 2) (c e k 3 ,C e k 4 ,C a ki,C a k 2 ,C a k 34 ) 



T T r 3 r 3 



Z b (c a ki,Cak 2 ,C e k 3 ,C e k 4 ,C a k 12 ) + 22^ (c a ki,C a k2,C e k 3 ,C e k 4: ,Caki2) 



(3) 



2T^c*CgI 6 (3) (c a ki,c e k 3 ,c a k 2 ,c e k 4 ,c e ki 3 ) 



2T s T c c 3 c 3 



Z b (c a /c2,c e fc 4 ,c a fci,c e A; 3 ,c e fe24) +Z b (c a ki,c e k 37 c e k 4 ,c a k 2 ,c e k 13 ) 



-T 2 c 4 c 2 

2 c^e^a 



Zfc 1 - 1 (c Q fci,c e fc 3 ,c a fc 2 ,c e fc4,Cefci2) + 2J^ 2) (c a fci , c e fc 3 , c e fc 4 , c a k 2 , C e fcl 3 ) + Z^ (c e k 3 , c a k 1: c e k 4 , c a k 2 , c e k 13 ) 



- 3H CT T s c^Cgl| 1 - ) (c Q fci,c a fc 2 ,c e A: 3 ,c e fc4,c a fci2) - 3S CT T c c^CeX( 2) (c a fci, c Q fe 2 , c e k 3 , c e ki, c a k 12 ) 

- Z C Y < T cjcll ( c 1) (c e k3,c e ki,c a ki,c a k 2 ,c a ku) - 2H c T CT Cg4x| 2) (c e fe3, c e fc 4 , c a ki, c a k 2 , c a ku) 

- 4S c T 5 .CgC^ 2) {c a k u c e k 3 ,c a k 2 ,c e k il c e ki 3 ) 

- 2~ c ? c cl4 Z« ( k 2 ,c e k^c e k 13 )+l<p {c a k llCe k 3 ,c e k A , c a k 2l c e fci 3 )J + 23perms| , 

(3.6) 

where in the first line "5 perms" denotes the other 5 possibilities of choosing two momenta for Q a and two momenta for 
Q s out of the four external momenta. Note that in the permutations, the speeds of sound c a and c e are always associated 
with the given extra momenta. The purely entropic four-point function (Q 4 ) is 



(Q s (t, fei) Q 8 (t, k 2 ) Q s (r, fc 3 ) Q s (r, fc 4 )) 



SE 



4 1 1 

--(2ir) 3 5 3 (^2 k i)^^ 2 c c 8 e c 2 a 2 a (ceki, c e k 2 , c e k 3 , c e fe 4 , c a k 12 ) + -T^CgC^ 1 ' (c e fei, c e fe 2 , c e k 3l c e k 4 , c a k 12 ) 



+ T s T c CgC Q X^ 2) (c e ki,c e k 2 ,c e ks, c e k 4 ,c a k 12 ) + -T 2 c c 6 e lj; 3) (c e ki,c e k 2) c e k 3 ,c e k 4 , c a k 12 ) 



- E c T a c%c%XP (c e fci, c e k 2 , c e k 3l c e k 4l c a k 12 ) - Z c Y c c 7 e c a l( 2) (c e k 1 ,c e k 2 ,c e k 3l c e k 4l c a k 12 ) + 23perms| 



(3.7) 



3.2 Full Trispectrum for the Curvature Perturbation 

The scalar field perturbations Q a and Q s themselves are not directly observable. What we are eventually interested in is 
the curvature perturbation 1Z. As has been investigated in $5% |38| , p3| l, the comoving curvature perturbation 1Z is related 
to the adiabatic and entropy perturbations of the scalar fields by 







I . 1 Qcr* 









(3.8) 



= NaQcr* +NsQs 



Here T-jis is the so-called transfer function from entropy perturbation to adiabatic perturbation 8 . Note that Af s = T-jisJ\f a 
is in general time-dependent. Thus contributions to the four-point correlation function for 7Z are given by 



{1l{k 1 )'R{k 3 )'R{k 2 )'R{k 4 )) 
■ A/" 4 [(Qa(ki)Qa(k 2 )Qa(k 3 )Q a (k A )) SE + (Q CT (fci)Q a {k 2 )Q a {k 3 )Q a (k 4 )) c ] 
NlN 2 [(Qa{ki)Qa(k 2 )Q s {k 3 )Qs(k 4 )) SE + {Qa{ki)Qa{k 2 )Q s {k 3 )Qs{k 4 )) c + 5 perms] 
+ Ns [(Q s (ki)Qs(k 2 )Qs(k 3 )Qs(k 4 )) SE + (Q s (fei)g s (fc 2 )Q s (fc 3 )g s (fc 4 )) c ] , 



(3.9) 



where the subscript g e denote the four-point functions which come from exchanging intermediate scalar modes, and the 
subscript c denotes the four-point functions which come from contact diagrams. The four-point functions of exchanging 



8 As was pointed in Jl la, |l 17j, as long as the fields roll slowly, these additional contributions after horizon-crossing are heavily suppressed. 
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scalar modes for Q a and Q s are given by (3.5), (3.6) and (3.7). The four-point function of the contact diagram is given in 



Ref.[104]. In deriving (3.S), we have used the assumption that there is no cross-correlation between adiabatic and entropy 
modes, i.e. (Q a Q s )* = 0, around horizon-crossing. 
It is convenient to define a so-called trispectrum 

4 



i=l 
4 

V '■ X " 

i=l 



(2n) 3 6 3 ( V ki){T c (kx,k 2 , k 3 , k 4 ) + T s (k u k 2 ,k 3 , fc 4 )) 



(3.10) 



where 7~ c is given by eq.(5.29) of Ref.[104|. We can derive T s from eqs.(3.5)(3.6)(3.7). 



To investigate the size of Non-Gaussianity roughly, we choose regular tetrahedron limit, k\ = k 2 — k 3 = k 4 = k\ 2 = 
k 3 4, and take the approximation c 2 — c 2 = c 2 <C 1. We define a real number t^L from the trispectrum to characterize its 
size, 



T(ki,k 2 ,k 3 ,k 4 )\ rt h = PntNL 



We have 



(3.11) 



tNL = t C NL +t 



where t c NL comes from the contact diagram [104] 



(3.12) 



*NL = (1 + Tig) 3 (tl +t 2 + t 3 ) 



where 



3c 2 (54c 2 A 2 - H 2 e(3\ + 10U)) 
512H 4 e 2 



-T4 



9c s (H 2 e - 15cf A) 



r2, i q„4> 



*3 



256c5i/ 2 e 
13 



T, 



81 



ns 1024c? 



l 3 (_ g 2 £ + 3c 4 A ) 13(-g^ + 3efA) ■ 

256c?tf 2 e RS 128dH 2 e ns 256c! 



515 

8192c? 



103 

2048c? 



comes from scalar exchange diagram, 



1-nl — + ts + te , 



where 



h 

*6 



2ct 



H 2 ( 



0.22c: 



0.67A\ 



AS 



0.06c7 4 T 4 



as > 



2.74c7 4 + (8.53 + 12.95£ + 5.75£ 2 ) c~ A T? 



RS 



-2.25-^- + (20.61c; 4 + 14.72c7 4 £ + 2.21^ 



1.43 + 1.99£ + 1.24£ 2 ) c~*T RS 
A 



2.30C 



H 2 t 



(3.13) 



(3.14) 



(3.15) 



Tj^ + (0.37 + 0.38£)c; 4 T^ s , (3.16) 



Here the contributions t±, t§, t§ come from diagram I a , and I c respectively (see Appendix |B| for details). Comparing 
ti, t^, te with ti, t 2 , t 3 , we can see that the scalar exchange diagram makes a nontrivial contribution to the trispectrum. 
As for the contact contributions, the contributions to the trispectrum from exchanging scalar modes can be enhanced by 
small sound speed(s), large Tus, large £, and large jh-. 

4. Conclusion 



In this note, we made a complementary calculation of the contributions to the trispectrum of primordial curvature pertur- 
bations from exchanging intermediate scalar modes in the context of generalized multi-field inflation, which completes 



the calculation of our previous investigation [104]. We choose regular tetrahedron limit to estimate the size of non- 



Gaussianity. The calculation presented in this work, together with [ 104 ], can be employed as the starting point for further 



analysis of the trispectrum of generalized multi-field inflatioyn models, such as the shapes, squeezed limit [ 1 09 , 110, 111] 
and estimators [ 1 12| , 113 , 1 14 , 115 , 118 ] etc. We would like to come back to these issues in the near future. 
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A. Coefficients in the interactional Hamiltonian 



The variously introduced coefficients in (3.3 ) are given by 



4A 



T, 



_jWe_ 

V2XP. 
1 



x \cl-l 



1 



XP 



x 



Hs/Te W a - 1 



V2 ( 



1 



(A.l) 



B. Basic Integrals 



The full expressions for the four-point functions are rather complicated. In this work, at the leading-order, all contributions 
to the four-point functions can be grouped into six basic inte 
and their "conjugate" which we define as below (see Fig. |b]). 



to the four-point functions can be grouped into six basic integrals, which we denote as I a , 7, (1) , 7, (2) , 7 fc (3) , iP and 7 C (1) 




r(i) 



r(2) 



r(3) 




la 



Figure 1: Diagrammatic representation of the six basic integrals: I a , J 6 (1) , 7< 2) , if \ and iP 



All 



the momenta configurations and n, T2 are the same as in I a - A red dot denotes the temporal derivative, a 
blue dot denote the spatial derivative or momentum in Fourier space, where a blue line between two dots 
represents the dot product. 



From Fig.^, it is straightforward to read the expressions for these integrals, we find 



I a {k 1 ,k 2 ,k 3 ,k 4 ,k 5 ) 
1 ft 



2H 2 
1 

4ff2 



/ dn / dT 2 -^—diG kl (r, n) diG k2 (r, Tx)d 2 Gk 3 (r, r 2 ) d 2 G ki (r, t 2 ) di 2 G k& (n, r 2 ) 

J-oo J-oo r l T 2 

\-r [ dn [ dr 2 -^—diG kl (n, r) diG k . 2 (n, r) d 2 G k3 (r, r 2 ) 9 2 Ga; 4 (r, r 2 ) <9i 2 G fc5 (n, r 2 ) , 

H J-oc, J-oc T \T 2 



(B.l) 
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where and in what follows d\ t2 = ^777^ $12 = dT f dT2 an d i n this appendix we denote G k (n, r 2 ) = u k {rn)u* k {ri 2 ) with 



If' (k 1 ,k 2 ,k 3 ,k4,k 5 ) 

(fci • hi) (fe 3 • M k 



2ff 2 



+ Iff 2 ( fcl ' ^ (** ' * 4 ) 
7 b (2) (fci,fe2,fc3,fc4,fc 5 ) 

(fci • fe 2 ) (fe 5 • fe 4 ) 



2i7 2 



/ (in / dr 2 — !— G kl (r, n) Gfc 2 (r, n) G fc3 (r, r 2 ) G fe4 (t, ^2) <9i 2 G fe5 (n, r 2 ) 

J-00 J-00 T 1 T 2 

r 1 

rfn / dr 2 G kl (ri,r) Gfc 2 (n, r) G fe3 (r, r 2 ) G fe4 (r, r 2 ) <9i 2 G fc5 (n, r 2 ) , 

X) J-CO T lT 2 

(B.2) 



rfn / dr 2 Gfe x (t, n) Gfe 2 (r, n)d 2 Gk 3 (r, r 2 ) Gt 4 (r, r 2 ) <9iG/c s (n, r 2 ) 

, TiT 2 



+ 71H ' fe 2) (fes • k 4 ) [ dn [ dr 2 -^—G kl (ti,t) G fe2 (n, r) d 2 G ks (r, r 2 ) G fe4 (r, r 2 ) <9iG fe5 (n, r 2 ) , 

4 # J-co J-co T 1 T 2 



J & (3) (fci,fe 2 , fc 3 ,fe 4 ,fe5) 



(B.3) 



1 



~2# 2 
~ Iff 2 



(fe 5 • fc 2 ) (fe 5 • fc 4 ) sft 



/ dn [ dT 2 ^—dxG kl (t, n) G fc2 (t, n) <9 2 G fc3 (r, r 2 ) G ki (r, r 2 ) G fcs (n, r 2 ) 
.J-co J-co nr 2 

^-r (fe 5 • k 2 ) (k 5 ■ k 4 ) I dn [ dr 2 -^—diG kl (ri,r) G k2 (n,r) <9 2 G fe3 (r, r 2 ) Gfc 4 (r, r 2 ) G fe5 (ri,r 2 ) . 

H J -co J— 00 r ir 2 



(B.4) 



and 



4 1} {ki,k 2 ,k 3 ,k 4 ,k 5 ) 

1 /" T /" Tl 1 

= 7772 ( fc 3 ' fe 4) 5ft / dn / dr 2 diG kl (r, n) <9iG fc2 (r, n) G fe3 (r, r 2 ) G fe4 (r, r 2 ) d 12 G ks (n, r 2 ) , R „ 

2-n J-co J-00 T ir 2 (.d.j; 

~ 777? ( fe 3 • fe4 ) / dTl / dT 2 -^—d x G kl (n, r) aiG fc2 (n,r) G fc3 (r, r 2 ) G fe4 (r, r 2 ) <9i 2 G fc5 (n, r 2 ) , 

4 -" 7-co J-co TlT 2 

7 C (2) {k 1 ,k 2 ,k 3 ,k 4 ,k 5 ) 
1 



= 2ff2 ^ 5 ' fc4 ) ^ y dri J dT2 ^T dlGk i ( T ' T 9iGfe 2 (t, ti) d 2 G k3 (r, r 2 ) G fe4 (r, r 2 ) <9iG fe5 (n, r 2 ) ^ ^ 
- 777? ( fe 5 ' fc 4) / (iri / dT 2 -^—dxG kl (n, r) <9iG fc2 (n, r) <9 2 G fc ., (r, r 2 ) G fe4 (r, r 2 ) <9iG fe , (n, r 2 ) . 

J_ !x> J_ ca TlT 2 

It is useful to introduce the "conjugate" contributions, defined as follows. Up to the second-order in perturbation 
theory, there are two interaction vertices and thus two temporal integrals with respect to n and r 2 respectively. We 
call two contributions (diagrams) are conjugate to each other with exchanging n «-> r 2 while keeping all the momenta 
relations. Having known the expression for a diagram, it is easy to write down the integral expression for its conjugate, 
e.g. 



I a (ki,k 2 ,k 3 ,k 4 ,ki 2 
1 



2H 2 
1 

Iff 2 " 



If 



/ dn I dr 2 -^—d 2 G kl (r, r 2 ) d 2 G k2 (r, r 2 ) <9iGa; 3 (r, n) diG ki (r, n) di 2 G k& (n,T 2 ) 

J— 00 J — 



00 r l T 2 
1 



(B.7) 



rfr 2 

co r l T 2 



9iG fel (n, t) <9iG fc2 (n, t) 9 2 G fe3 (r, r 2 ) 9 2 G fc4 (r, r 2 ) 9i 2 G fe5 (n, r 2 ) 



where * denotes complex conjugate. It is analogous for the other conjugate integrals, which we do not write here for 
simplicity. Moreover, we introduce the combination of a contribution and its conjugate, e.g. 



l a {ki,k 2 , k 3 , k 4l k 5 ) = [I a + i a ] {ki,k 2 , k 3 , k 4l k 5 ) . 



(B.8) 



Before we evaluate the integrals, it is useful to make it clear about the smallest set of integrals we need. There are two 

(2) ~(2) (2) (2) "(2) 

cases. For left-right asymmetric diagrams, e.g. 1^ (or I '), we always encounter the combination!^ = I ' + 1^ 
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~<2) 

rather than If, itself. While for the left-right symmetric diagrams, e.g. I a , I a is simply exchanging simultaneously 
fci o fc 3 , fc 2 &4- Thus, after the 6 permutations (which specify two momenta associated with n and other two 
momenta associated with r 2 ) among the four extra momenta fci , • • • , k&, the final contribution to the correlation function 
from I a is equal to X tt /2. Thus, what we really need is the following six basic integrals: X a , I^K \ 1^ and 

T (2) 

Now we collect the final results for these integrals, in the limit of t — > 0. We find 



I a (k 1 ,k 2 ,k 3 ,k4,k 5 ) 



H s k 5 
1QK 5 



rH 



A (^12, if 3 4,fc 5 ) + 



K 5 



(#12 + fc 5 ) (#34 + fc 5 r 



(B.9) 



with 



_ 10g| + ( Sl + 3r) {As 2 + Jf) + 6r 2 
A(si,s 2 ,r) = -3 h (si <h> s 2 ) , 



(S2 + r-) 

where here and in what follows we denote = fcj + fcj and K = ki + k 2 + k 3 + fc 4 . 



(B.10) 



I b (1) (fci, fe 2 , fe 3 , ki, k 5 ) = (fci • fc 2 ) (fc 3 • fc 4 ) 



64if 5 



ni 



[r (#12, ^34; fcifei, fc 3 fc 4 ; J12) + (12 34) + tf 5 F (ifi 2 , k 5 ,hk 2 ) F (K34, k 5 ,k 3 k 4 )] . 

(B.ll) 



with Jij = Kij — k 5 , and 



F(s 1 ,s 2 ;q 1 ,q 2 ,t) 



{K 6 + K 5 {-2t + Sl + 2s 2 ) 



\K-tf 

+ K 4 [t (t - 2si) + 3 (-t + si) s 2 + 2gi + 6q 2 ] 
+ K 3 [t 2 (si + s 2 ) + 8s 29l + 12si 92 - t (5sis 2 + 4 ft + 6q 2 )] 
+ K 2 [2t (ts lS2 + (t - 7s 2 ) qi ) + 2(t 2 - 8isi + 20 qi ) q 2 ] 
+QKt [ts 2 qi + {tsi - 10<7i) q 2 ] + 24t 2 q 1 q 2 } , 



(B.12) 



and 



F(s,t,q) 



2s 2 + 2q + 3st + t 2 
(s + tf 



(B.13) 



l£ 2) (fci,fc 2 ,fc 3 ,fc 4 ,fc 5 ) 



= (fci • fc 2 ) (fc 4 • fc 5 ) — 



h\ 2 H s k 5 (X 1 



k 5 J UK 5 



n fc 3 



\i=l 



x [r (i^i2, if 45 ; kik 2l k 4 k 5 ; J i2 ) + T (if 45 , K 12 ; -/c 4 /c 5 , fcife; J 34 ) + K 5 'F (K 12 , k 5 , k x k 2 ) F (K 45 , k 3l fc 4 fc 5 )] , 

(B.14) 

where = ki — kj . 



Xt 3) (fci,fc 2 ,fc 3 ,fc 4 ,fc 5 ) 



: (-fc 2 • fc 5 ) (fc 4 • fc 5 ) 



k-y k^ H 
"kfUK 5 



kf 



x [r (K 25 ,K i5 ; -k 2 k 5 ,k±k 5 ; J 12 ) + T (K 45 , K 25 ; -fc 4 fc 5 , fc 2 fc 5 ; J 34 ) + K 5 F (K 25 ,ki, k 2 k 5 ) F (K 45 , k 3 , fc 4 fc 5 )] . 

(B.15) 



And 



( k i,k 2 ,k 3 ,k 4 ,k 5 ) 

E(fc 3 -fc 4 ) 



32^ 5 k x k 2 klk\ 



C(K 34 , k 3 k 4 , J12) + C{K 3il k 3 k i7 J 34 ) + 



if 5 F(g 34 ,fc 5 ,fc 3 fc 4 ) 
(^i2 + fc 5 ) 3 



(B.16) 
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1 {2) {k 1 ,k2,k 3 ,k 4 ,k 5 ) 

(B 



H 8 1 

= (fc 5 • fe 4 ) 



32K 5 k^hklh 

with 



C (K 45 , k,k 5 , J 12 ) + C {K 45 , -hk 5 , J 34 ) + K F {Ki5 ' **' hh 



(k 12 + h) 

K(K-t) [-t (K + 3s) + K{K + 4s)] + 2 (lOK 2 - 15Kt + 6t 2 ) q 



C(s, q ,t)-— ^-^ - 

^(3^ (if + 2s) +t 2 (if + 3s) - ift (3/f + 8s)) +2 (ldfiT 2 - 15ift + 6i 2 ) g <B 
C (s, q, t) = — ^ — . 

(K-tf 
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